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ABSTRACT 

Aims. We investigate the anisotropy of stellar winds in binaries to improve the models of accretion in high-mass X-ray binaries. 
Methods. We model numerically the stellar wind from a supergiant component of a binary in radial and three-dimensional radiation 
hydrodynamic approximation taking into account the Roche potential, Coriolis force, and radiative pressure in the continuum and 
spectral lines. 

Results. The Coriolis force influences substantially the mass loss and thus also the accretion rate. The focusing of the stellar wind by 
the gravitational field of the compact companion leads to the formation of a gaseous tail behind the companion. 
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1. Introduction 

Roche-lobe overflow and an isotropic radial stellar wind have 
traditionally been considered as two competing scenarios to ex- 
plain the mass loss from a companion star that can feed an ac- 
cretion disc around a compact object, particularly in interact- 
ing b inaries and X-ray binaries - cf., e.g. Shaku ra~& Sunyaev] 
( |1973| ). One of us thus developed a model of an evaporative ra- 
dial stellar wind modulated by the effective Roche potential to 
overcome this dichotomy (Hadra va|1980| ). This model revealed 
that there is a smooth transition between the outflow from a nar- 
row stream around the inner Lagrangian point L\ for small ra- 
tios of thermal to escape velocities, to an almost isotropic wind 
for high values of this ratio. The distribution of mass-loss rates 
on the surface of the companion is modulated depending also 
on the mass ratio of the system. In some cases, the mass loss 
in the vicinity of the point exceeds that around L\. S. Kffz 
(1980, private communication) raised an objection to this model 
based on the argument that it neglects the Coriolis force and 
thus artificially enhances the outflow by an unrealistic centrifu- 
gal force. Because of this criticism, the model was not published 
until Hadrava (1987 , cf. the on-line Appendix). In the meantime, 
a similar model of a r adiatively driven wind was published by 
|Friend & Castor| ( [l982| ). This model of a so-called focused stel- 
lar wind is now widely used to explain the behaviour of X-ray 
binaries. However, the above-mentioned objection is also valid 
for this model and the Coriolis force should be taken into ac- 
count to improve the model of stellar wind anisotropy. This may 
significantly influence the interplay between the mass loss from 
the donor star, its accretion by the compact component, and 
the consequent ionization of the dono r's wind and atm osphere 
illuminated by the X-ray radiation - cf. |Gies et al.| ( |2008] ). 

Here we present our first preliminary results of our revised 
calc ulations of the anisotropic stellar wind in binaries achieved 
by |Cechura| ( |20 1 0| . In Section|2j we summarize our calculations 
of the radial wind based on the Friend-Castor approximation. 
To reveal the role of the Coriolis force, we then numerically 



modelled on a three-dimensional (3D) grid the radiation hydro- 
dynamics of the stellar wind taking into account the same ex- 
pressions for the radiative pressure in the lines - cf . Section [3] 
These calculations showed the formation of a gaseous tail be- 
hind the compact companion, which resembles Bondi-Hoyle- 
Lyttleton accretion (Hoyle & Lyttleton 1939 ) - cf. Section [4] 



As discussed in Section |5J this feature actually deserves to 
be denoted as a gravitational focusing of the wind rather than 
the slight tidal modulation of the wind's root. 

2. Radial stellar wind 

Our present model of anisotropic radial stellar wind is based 
on a modified version of the line-driven wind theory of Castor, 
Abbott, and Klein ( [Castor et aL]|1975| CAK hereafter), where 
the wind material is predominantly accelerated by the line ab- 
sorption of the star's radiation field. In this work, the authors 
solved the equation of radiative transfer in the stellar wind under 
the assumption of the Sobolev approximation ( Castor|1974 ) and 



found that the line force could be approximated by a power law 
in optical depth 

where t is the optical depth parameter defined by 

dvY 1 



t = aQpVth [drj ' (2) 

and k and a are parameters of the CAK model that depend on 
the effective temperature of the star. Physically speaking, k is 
a measure of the number of driving lines, and a is a measure 
of the relative distribution of optically thick to optically thin 
lines. We chose the values of a = 0.7 and k = 1/30 in our 
model. These values should provide a reasonable approximation 
of the line force corresponding to a star with an effective temper- 
ature of around 30000 K (cf. CAK). The quantity cr e is the elec- 
tron scattering coefficient, L* is the luminosity of the primary 



1 



P. Hadrava and J. Cechura: Really focused stellar winds in X-ray binaries 



star, p is the density of the wind material, and v± is the thermal 
velocity of the ions in the wind. 

In the case of an isothermal radial outflow, the equation of 
momentum conservation is given by 



.x10 
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(3) 



where a = (dP g /dp^ is the isothermal sound speed, M* is 
the mass of the primary star, and T* is the ratio of the primary 
star luminosity to the Eddington luminosity 



r* = 



4ttGM*c 



(4) 



To include the effects of the compact companion on the wind 
dynamics, we have replaced the spherically symmetric gravita- 
tional potential of the mass-losing star by the Roche potential 
® e ff , i.e. both the tidal force and the centrifugal force due to a ro- 
tation synchronized with the circular orbital motion of the com- 
panion are taken into account. In a spherical coordinate system, 
a new effective potential is given by 



O cff (r,0,0)= - 



GM (D(l -r*) 
D 

-r x ) 



[l-2(r/D)A + (r/D) 2 ] 
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where the variables A and ji are defined by 

A = cos cp cos 6 , 
\i - sin0 , 



(5) 



(6) 



in terms of the latitude and the longitude 6 on the surface 
of the mass-losing star, D is the separation of the components 
of the binary, q is the ratio M x /M*, where M x is the mass of 
the compact companion, and = and 6 = determines 
the line-of-centres of the two stars in the orbital plane. 

The radiation pressure in the continuum caused by 
the Thomson scattering of the photons emitted by both com- 
ponent stars on free electrons is also involved by means of 
the Eddington factors T* and T x , which are chosen as free pa- 
rameters for both stars (for which we use the values 0.26 and 
0.31, respectively). 

These forces cause a shift in the critical point of the outflow 
that varies across the surface of the mass-losing star and they 
distort the spherical symmetry of the wind. The mass-loss rate 
becomes a strong function of the size of the primary component 
relative to its critical potential lobe. In this radial approxima- 
tion, we still ignore the tangential forces caused by the pres- 
sure gradients between the neighbouring flux-tubes, as well as 
the Coriolis force, which should be taken into account when 
computing in the non-inertial reference frame co-rotating with 
both stars. The streamlines of the wind are strictly radial and 
the material is confined within a selected cone, hence the mass 
continuity equation becomes 



dM 
dQ 




■ pvr 



(7) 
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Fig. 1. Distribution of stellar wind intensity on the surface of 
the optical component of Cyg X-l in the radial approxima- 
tion (CAK model - the uppermost panel), the radial approxima- 
tion of the radiation-hydrodynamic model in 3D Cartesian grid 
(the second panel), the non-radial radiation-hydrodynamic so- 
lution without the Coriolis force (the third panel) and the full 
radiation-hydrodynamic solution with the Coriolis force in- 
cluded (the bottom panel). The direction to the L\ -point and 
the companion has the sidereographic longitude 0°, the longitude 
180° is directed towards the L2 -point. The number in the lower 
right corner of each panel gives the overall outflow of the mate- 
rial from the supergiant in M Q yr~ l 



where dM/dQ is the mass-loss rate per unit solid angle. The out- 
flow of the wind is artificially reinforced by an unrealistic cen- 
trifugal force. We expect this approximation to have little influ- 
ence on the deep highly subsonic layers of the wind, although 
it must lead to a completely unrealistic model in the supersonic 
circumstellar region. 

The results of this radial numerical model are presented in 
the upper panel of Figure[T] where we display the directional de- 
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pendence of the stellar wind intensity (i.e. the mass-loss rate per 
unit solid angle). We have chosen values of these parameters that 
are appropriate for the high-mass X-ray binary system Cygnus 
X-l and its primary O star supergiant HDE 226868. The orbital 
period is 5.6 days and the masses of the supergiant and the black 
hole have been chosen to be equal to 24. and 8.7 M©, resp., 
within the still wide region of the observational limits 23*? M© 



Caballero-Nieves et al. 



(2009). These val- 



and ll^M© found by 
ues result in a distance ot 42.4 K Q between the component cen- 
tres. The critical Roche lobe of this system has a mean radius 
of 20 R Q : the L\ -point is at a distance of 25.6 7?©, while the in- 
tersection of the critical equipotential with the rotational axis 
is 18.7 R Q from the supergiant centre. We have set the surface 
of the donor star to a mean radius 18 7?© (i.e. at the equipoten- 
tial of size 19.9 R Q towards the companion and 17.3 R Q towards 
the pole). The gravity acceleration g varies across this equipoten- 
tial from 1 .5 to 3. 1 x 10 3 cm.s~ 2 towards the L\ and the pole, resp., 
i.e. the mean value lo g g = 3.3 is close to the o bservational limit 
3.00 ± 0.25 found by |Caballero-Nieves et al.| p009). Similarly, 
the temperature 30.0 kK, which we have taken from the values 
for which CAK gives the parameters k and a in Eq. ([T]) , roughly 
corresponds to the value 28.0+2.5 kK determined in Caba llero- 
INieves et al.] ( |2Q09] ). 

The proximity of the black-hole companion and the assumed 
synchronous rotation of the mass-losing supergiant obviously 
disturb the spherical symmetry of the wind and the most in- 
tensive outflow of the material is concentrated in the equato- 
rial plane. The highest intensity can be found at a point directly 
facing the black hole. However, more interestingly, another sig- 
nificant local maximum is situated on the completely opposite 
side of the star. According to our model, the total mass-loss 
rate M of HDE 226868 is close to 2.23 x lO^M©^" 1 . The to- 
tal outflows from the facing and opposite hemispheres are al- 
most comparable and their ratio depends on the system parame- 
ters. These results qualitatively agree with the earlier works of 
|Friend & Castorl ffe82] ) and |Hadrava| ( [T987] ). It should, how- 
ever, be noted that there is a significant quantitative difference 
between the tidal modulation of the evaporative and line-driven 
winds. The typical thermal velocity is approximately 22 km/s for 
the temperature 30 kK, which, compared to the critical rotational 
velocity 500 km/s, gives a ratio of the thermal to the binding 
potential energy of the order of 10" 3 . The Parker's evaporative 
wind is thus highly concentrated in the vicinity of the L\ -point 
(cf. the Appendix |A.2| Fig. |A.4| ). However, if the CAK- radiative 
drag enhances the wind and the radiation of the star is assumed to 
be isotropic, the tidal modulation of the wind intensity is partly 
smeared out. The distribution of mass-loss rate is then sensitive 
to the choice of the surface level. 



3. Radiation hydrodynamic model of the wind 

3.1. Equations of motion 

To deal with the previously mentioned objection to the radial 
approximation, we also created a numerical model of the stellar 
wind based on the radiation hydrodynamics calculated in a 3D 
Eulerian coordinate grid. In this case, a new equation of motion 
is given by 

5 = -(" v > " v °eff + Sl ~ ~ VP g + 2v An , (8) 
ot P 

which has to be solved together with the continuity equation 

dp 



at 



+ V • (pv) = o . 



(9) 



The constant orbital angular velocity n is given by the third 
Kepler's law (n 2 D 3 = (1 + q)GM*). Our present model uses 
the purely Newtonian approach to handle the gravity and the dy- 
namics. We can naturally expect that relativistic effects are 
important in the close vicinity of the compact companion in 
the X-ray binary, which is either a neutron star or a black hole. 
However, the grid we use for our numerical solution is sparse 
in order to achieve a reasonable computational time. A typi- 
cal dimension of one cell corresponds to several thousands of 
Schwarzschield radii and we are unable to distinguish structures 
that are smaller than this typical length. It is therefore reason- 
able to assume that at such a distance, relativistic effects will be 
negligible. 

We have assumed that the radiative pressure /l in the lines, 
which is a main source of the repulsive force, is radial and that 
its value is also given by Eq. ([I]). The gravity and continuum 
radiation pressure of both components of the binary have also 
been included in terms of the effective potential ® e ff given by 
Eq. ([5]). In addition, new effects are taken into account, such 
as the tangential components of the gas-pressure gradient and 
the Coriolis force. The gas pressure P g was taken from the equa- 
tion of state of an ideal gas whose temperature is fixed to the ef- 
fective temperature of the star in the present version of our code. 
This is a simplifying approximation that we will improve in fu- 
ture versions of the model, although it may not be too inaccu- 
rate in the case of a rarefied gas in thermodynamic equilibrium 
with the radiation of both stars. To develop a self-consistent 
radiation-hydrodynamical model, the radiative transfer should 
also be solved, whereas in our present model the radiation field 
is prescribed. In the subsequent paragraph, we briefly discuss 
the possible consequences of this simplification. 

3.2. The radiation field 

In his referee report to this paper K. Gayley raised additional 
objections to the Friend-Castor approximation used to describe 
the radial line-driven wind in binaries. On the basis of their in- 
vestigations of the stellar winds of early-type fast rotators and 
the problem of the (non-)existence of wind-compressed disks 
( |Owocki et al.| 1 996 1 [T998| |Gayley & Owocki|2000| ) Gayley sug- 
gested that both the non-radial components of the radiative force 
and the anizotropy of the supergiant radiation due to the gravity 
darkening may substantialy influence the structure of the wind 
in binaries. To verify whether this suggestion is indeed valid, 
we would need to create a sophisticated model that exceeds 
the possibilities of our current numerical model and intentions 
of this work. We thus give here only our preliminary estimates 
of the possible influence of these effects. 

The non-radial components of the radiative force in super- 
sonic axisymmetric flows were studied by Grinin (1978). In a 
differentially rotating medium such as a Keplerian disk, they 
are due to the asymmetry of the velocity gradients in directions 
skewed forward and backward with respect to the radial direc- 
tion and the consequent difference in the depletion or enhance- 
ment of the line radiation intensities. As for either the Poynting- 
Robertson effect ( Robertson|1937|) or the reflection of radiation 
by rotating disks (Hadrava et al. 1997), this radiative torque can 
transfer the energy and angular momentum from the medium to 
the radiation and stimulate a fall down of the flow. Gayl ey &] 
Owocki (2000) found that this effect may decrease the angular- 
momentum loss of the fast rotators by the wind to about 30% or 
40%. 

However, in the case of the synchronously rotating compo- 
nent of a binary with a non-negligible mass of the companion, 
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the rotation is safely subcritical (even if the component fills its 
Roche lobe, which is restricted by the tides). In the particular 
case of our present model, the rotational velocity of the super- 
giant is approximately 160 km/s, while the critical velocity is 
500 km/s. Consequently, the influence of the non-radial compo- 
nents of line-driven radiative force is also smaller than in the case 
of critically rotating stars. 

While these non-radial components of the radiative force 
may to some extent influence the motion of the circumstellar 
matter in the supersonic region, where a non-local coupling can 
also take place, their role is even more questionable in the sub- 
critical region where the surface distribution of the wind inten- 
sity is determined. In this quite narrow slab of the atmosphere, 
the tangential component of the wind velocity increases only 
slowly due to the Coriolis force, hence the optical depths in 
the skewed directions given by Eq. ^ cannot differ significantly. 

A more substantial objection to the Friend-Castor model is 
the possible r ole of g ravity darkening. According to von Zeipel's 
theorem (von Zeipel|1924 ), the distribution of the radiative flux 
F and the effective temperature across the tidally or rota- 
tionally distorted surface of a star should be given by the local 
gravity acceleration g as 



eff 



(10) 



where J3 = 0.25 in the standard von Zeipel's formula. Because 
the radiative force f L given by Eq. ([T]) is proportional to the flux, 
it may enhance the wind more in the polar region of the star, 
where g reaches its highest value, than in the equator and espe- 
cially in the line joining the component stars, where it has its 
smallest value. For the values of parameters chosen in our cal- 
culations, g is higher by 23% at the poles and lower by 15% and 
40% in the directions towards the L 2 and L\, resp., than its value 
in the perpendicular direction. This effect is thus the opposite of 
the direct modulation of the wind by the effective potential. In 
the region, where the distance from the star is comparable to its 
radius, this effect is additionally enhanced by the ellipticity of 
the star. 

To determine the outcome of these mutually competitive ef- 
fects, one would need to develop a reliable model of the gravity 
darkening. The above-mentioned original von Zeipel's formula 
is widely used owing to its simplicity, but it is neither physically 
self-consistent nor confirmed by observations, which mostly in- 
dicate that the actual gravity darkening is smaller, as expected 
from a more detailed theoretical treatment. |Lucy | ( [T967] ) reinves- 
tigated the problem of gravity darkening asuming a convective 
instead of a radiative equilibrium arriving to a generalized for- 
mula in Eq. ( [T0| ) in which the exponent ft = 0.08. The classical 
von Zeipel's formula is based on the solution of radiative trans- 
fer in the diffusion approximation only, assuming hydrostatic 
equilibrium and hence the homogeneity of the atmosphere on 
equipotential surfaces (cf. Appendix |A.1| ). The simplified treat- 
ment of radiative transfer neglects both the dilution of the radia- 
tive flux in the parts of sub-surface layers with higher curvature, 
which enhances the gravity darkening, and the tangential dif- 
fusion of the radiation, which decreases its effect, (cf. Hadrava 
|1992| ). However, even more importantly, the resulting gravity 
darkening violates the homogeneity and the hydrostatic equilib- 
rium, which are the starting assumptions of von Zeipel's theo- 
rem, and gives rise to meridional circulations. The Coriolis force 
acting on these flows may then also break the mirror symmetry 
of the advancing and recessing hemispheres of the star and hence 
the simple picture of the flux given by the local value of g (cf . 
Gay ley et al.|2007| ). Our detailed solution of radiative transfer in 
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Fig. 2. Distribution of the stellar wind intensity in dependence 
on gravity darkening. The radial approximation for the same pa- 
rameters as in Fig. [T] is calculated for radiation modulated ac- 
cording to Eq. ( [T0| ) and values of the parameter J3 from the top to 
the bottom of 0.0, 0.125, 0.1875, and 0.25 



yields j3 ^ 0.13 for small variations in g and a more complex 
dependence F = F( g) for its larger vari ations (cf . Hadrav a[l987| 
and Appendix A.l particularly Figure \AJ). Yet another value 
P 0.19 was obtained by |Claret| p012), who also gives refer- 
ences to several observational results confirming that the actual 
gravity darkening parameter has a scatter and is usually smaller 
than the von Zeipel's value of 0.25. 

On the referee's suggestion, we recalculated the outflow 
in the radial approximation described in Section |2j assuming 
the gravity darkening given by Eq. ( [TO} with several values of 
the parameter J3 (cf. Fig. [2]). For the medium value fi = 0.125, 
the tidally enhanced outflows in the directions of the L\- and 
L 2 -points can really be diminished due to the decrease of the ra- 
diative drag. For even higher values of ft, our calculations show 
a formation of local minima of the wind intensity in these direc- 
tions, although the local minima caused by the potential barrier 
in polar regions persist and the maximum of the wind intensity 
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appears at moderate latitudes. This is, however, an upper esti- 
mate based on a model with several inconsistencies. 

Because the starting assumption of hydrostatic equilibrium 
is violated from the very beginning in the case of strong stel- 
lar winds, we have not included the von Zeipel's formula for 
gravity darkening in our 3D models. To achieve a really self- 
consistent radiative-hydrodynamic solution, the 3D radiative 
transfer should be included in the procedure, which would be 
computationally much more demanding even in the simpler case 
of axially symmetric rotating stars. In the current approxima- 
tion, we include in the hydrodynamic modelling the radiative 
forces from the given purely radial and isotropic radiation of 
the components and postpone the treatment of radiative transfer 
in the circumstellar matter to a future study. 

3.3. Numerical model 

The results of the previously computed radial model, i.e. the spa- 
tial distribution of density p and the Cartesian components of 
the velocity vector v, were used as initial conditions for the cor- 
responding quantities in the new 3D hydrodynamic simulation 
of the time evolution of the wind and circumstellar matter. 
They were also used to define the inner boundary conditions for 
the gas outflow in the photosphere of the donor star, which is 
nearly in hydrostatic equilibrium and the effects of the additional 
terms could thus be neglected there. 

An equidistant Cartesian coordinate grid was chosen in 
the co-rotating reference frame. The values of p and v found 
from the radial-wind model were interpolated onto this grid and 
kept fixed inside a Roche equipotential chosen to correspond to 
the surface of the donor star in order to define the inflow bound- 
ary condition for the Eulerian radiation-hydrodynamic solution. 
The outer boundary condition at the edge of the grid is a sim- 
ple outflow boundary condition with derivatives of all variables 
set equal to 0. The time derivatives of the density and compo- 
nents of the velocity vector were calculated from the equations 
of the mass and momentum conservation, resp. The small cor- 
rections corresponding to a sufficiently short time-step At were 
then added to the initial values. Repeating this process, the ini- 
tial stationary solution of the radial model adapted to satisfy 
the conditions of the new model. A new stationary solution was 
achieved after a relatively quick convergence of all quantities. 
The sum of time steps typically corresponds to a small number 
of orbital periods of the system (cf. Figure [3}. This is obviously 
a consequence of the wind velocity at the critical point given by 
the sound speed a being comparable to the orbital velocity nD 
in our investigated cases when the mass loss by stellar winds is 
substantial. The convergence could be much slower in the sub- 
sonic layers of the wind, if the inner boundary conditions were 
chosen on a deep equipotential where the density is high and 
outflow velocity low. However, to avoid the steep density gradi- 
ents in these layers, which could not be represented with suffi- 
cient accuracy in the limited space-resolution of the grid, we had 
to choose this layer closely below the surface of critical points. 
The convergence of the numerical model was thus accelerated by 
the numerical viscosity within a few grid cells around the inner 
boundary surface. 

A value of time step At, which is constant for the whole 
grid, was subjected to the Courant-Friedrichs-Lewy condition, 
specifically 

(1 Ax Ay Az\ 

-min(A*,Ay,Az),— , (11) 
a \v x \ \v y \ \v z \J 



where Ax, Ay, and Az are the distances between the neighbour- 
ing grid nodes in the x, y, and z-directions, resp. The Courant 
number Cq is a safety factor, which in this case is taken to be 
0.5. 

We chose a combination of the forward-time central-space 
scheme (FTCS hereafter) and the Lax-Friedrichs numerical in- 
tegration scheme in our computations. This combination allowed 
us to exploit all of the advantages of both schemes, in addition to 
suppressing most of their drawbacks. The FTCS is a very sim- 
ple and computationally fast integration scheme, which is un- 
fortunately always unstable. To overcome this instability, we in- 
troduced a correctional mechanism to prevent random numeri- 
cal peaks. On the other hand, the Lax-Friedrichs scheme is sta- 
ble in all cases when the Courant stability condition is satisfied 
and therefore it serves as a stabilizing component in our model. 
However, the Lax-Friedrichs scheme is characterized by a high 
numerical viscosity, which results in a high level of diffusion. 
For a large number of integration steps, this feature would lead 
to smearing of all steep gradients, even those that are physically 
interesting. This is why we applied the Lax-Friedrichs scheme as 
infrequently as possible, specifically only once in every 25 inte- 
gration steps. If we had included the Lax-Friedrichs scheme with 
a lower frequency, the numerical simulation would have become 
unstable in the long term. 

The additional averaging of the density and the velocity in 
the outer boundary of the integration area was added to prevent 
oscillations and numerical instabilities within the stationary so- 
lution. This averaging is responsible for less accurate results in 
this area. Nevertheless, the density is very low at the border of 
the integration area and therefore should not have a measurable 
impact on the solution in the inner regions. 

The evolution of the density distribution in the equatorial 
plane during the convergence process towards a new stationary 
solution is captured in Figure [3] The initial state corresponds 
to the conditions set up by the radial approximation. The inte- 
gration process was stopped when a nearly steady-state solution 
corresponding to the newly included forces - the Coriolis force 
and the tangential forces caused by the pressure gradients - was 
reached. The white index in the upper right corner of each panel 
indicates an evolutionary time expressed in the units of the or- 
bital periods. The red cross and the asterisk show the position of 
the black hole and the centre of the supergiant, resp. The patch 
at the centre of each panel represents the supergiant component 
of the binary. More specifically, it is an intersection of the vol- 
ume defined by an equipotential surface and the equatorial plane. 
The convergence to a new stationary solution is relatively quick: 
at the time index equal to 3.0, the changes of studied quantities 
in one time step are less than 10" 4 of their initial values. As indi- 
cated in Figure [3] the originally radial distribution of density be- 
comes strongly influenced by the presence of the compact com- 
panion. We can see a gradual formation of a density tail behind 
the compact component, which is diverted by the Coriolis force. 
It is also apparent that the density increases with time through- 
out the entire modeled space. The reason for this increase is that 
we reduced the centrifugal force by allowing the wind to be di- 
verted by the Coriolis force. Therefore, the repulsive net force is 
weakened, leading to the slowdown of the radial component of 
the wind velocity and hence due to the equation of mass conti- 
nuity also to the increase in the wind density. 

A modified directional distribution of the stellar- wind inten- 
sity is traced back to the roots of the streamlines on the Roche 
equipotential surface at the centre, which was selected to repre- 
sent a surface of the supergiant. The result for the same stellar 
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Fig. 3. Evolution of the density distribution in the equatorial 
plane computed with the radiation hydrodynamic model dur- 
ing the convergence process towards a new stationary solution. 
The uppermost left panel represents the state close to the con- 
ditions set up by the radial model. The bottom right panel cap- 
tures the nearly steady-state solution corresponding to the newly 
included forces - the Coriolis force and the tangential forces 
caused by the pressure gradients. The white index in the up- 
per right corner of each panel indicates an evolutionary time 
expressed in the units of the orbital periods. The red cross and 
asterisk show the position of the black hole and the centre of 
the supergiant, resp. 



parameters corresponding to the system Cyg X-l is shown in 
the bottom panel of Figure [T] 

To distinguish the role of the different numerics and 
the physical effects involved in the new model, we also cal- 
culated two artificial intermediate models: the radial model in 
the 3D Cartesian grid and a non-radial model with tangen- 
tial gradients and velocity components included but for which 
the Coriolis force had been neglected (cf. the second and 
the third panel of Figure [TJ resp.). 

The second panel of Figure [T] represents results of the 3D 
hydrodynamical simulation where the non-radial components of 
velocity were suppressed in each step of the computing process. 
We thus allowed the wind material to move in the radial direction 
only and forced it to reach a stationary solution independently in 
each radial cone. Consequently, this solution should correspond 
to the radial approximation depicted in the uppermost panel. 
However, the numerical results show a considerably higher in- 
tensity of the stellar wind all across the surface of the supergiant 
with a particularly enhanced outflow in the direction of the L\- 
point and the integrated outflow increase to 2.52 x 10" 6 M Q yr~ l . 



This difference between the two physically equivalent models is 
obviously due to numerical errors caused by the insufficient spa- 
tial resolution of the 3D-grid in the inner region of the comput- 
ing area where we meet a steep density gradient. As a result, we 
experience an increase in the density all across this region and 
consequently a higher level of intensity of the stellar wind. To 
achieve certain qualitative conclusions about the role of the non- 
radial motion of the wind and the Coriolis force in real binaries, 
we have to compare the subsequent results of our 3D hydrody- 
namical model with this normalized radial approximation solu- 
tion. 

In the third panel of Figure [T] the results of the midway hy- 
drodynamic simulation are presented. In this model, we also take 
into account the tangential forces caused by the pressure gradi- 
ents between the neighbouring flux-tubes. We set no additional 
conditions on the velocity so the streamlines of the wind can 
bend into non-radial directions though the Coriolis force is still 
excluded. We note that the overall intensity of the wind has in- 
creased compared to the normalized radial approximation solu- 
tion in the second panel. The directional distribution of the out- 
flow has also changed. The intensity in the equatorial plane has 
increased, in contrast to the drop-off that has occured in the po- 
lar regions. Allowing the non-equatorial streamlines to deviate 
from the radial direction, they can perceive the centrifugal force 
and are slightly bent towards the plane of the orbital motion of 
the binary. The most significant increase in the outflow at the pri- 
mary maximum happens because the streaming of gas through 
the Roche potential window around the L\ -point is more effec- 
tive if non-radial velocities are allowed. A larger amount of ma- 
terial is transported in the direction with the lower effective po- 
tential barrier and the overall intensity of the outflow increases 
to 2.72 x 10~ 6 M Q yr-\ 

3.4. Results of the computations 

The final solution of the 3D hydrodynamical model including 
the Coriolis force is shown in the bottom panel of Figure T] 
Comparing it with the results of the previously described model 
in the third panel we note a slight shift in the primary and sec- 
ondary maxima counter-clockwise in the equatorial plane. This 
effect can naturally be expected if we introduce the Coriolis 
force into our problem - most of the matter escapes along 
the streamlines that pass through the L\ or L2 points or their 
vicinities and originate in deeper layers shifted in the direction 
of the orbital rotation. Another important feature of the new sta- 
tionary solution is that the mass-loss rate at the primary maxi- 
mum is lower than the value of the same quantity in the previous 
model. We note that there is an even more significant decrease in 
the intensity at the secondary maximum on the opposite hemi- 
sphere from the black-hole companion. The overall intensity has 
decreased and the integrated mass-loss rate for the entire star 
has dropped to 2.50 x 10" 6 M Q yr~ l . This effect is also caused 
predominantly by the Coriolis force - the flux tubes originating 
in parts of the stellar surface where the gravity is most signifi- 
cantly reduced by the tidal force must overcome a higher barrier 
of the Roche potential, while the actual maximum of the outflow 
around the Lagrange points is less enhanced by the tides. In other 
words - we now allow streamlines of the stellar wind to curve 
and not to co-rotate with the orbital motion. In consequence, we 
reduce the centrifugal force, which repulses the material from 
the star, and therefore we may expect that the overall mass-loss 
rate be lower than in the hydrodynamical simulation in which 
the Coriolis force is not included, as predicted by S. Kffz. The ef- 
fects of the non-radial flowing and the Coriolis force practically 



6 



P. Hadrava and J. Cechura: Really focused stellar winds in X-ray binaries 




Fig. 4. Hydrodynamic solution of the wind-fed circumstellar 
matter. The white cross shows the position of the black-hole 
companion, the black crosses mark the positions of the points 
L\ and L^. The distance scale on the x-axis is labelled in 
units of the separation between both components of the binary. 
The right-hand-side colour-scale corresponds to the distribution 
of the wind density in the equatorial plane. The left-hand- side 
colour- scale describes the directional dependence of the inten- 
sity of the stellar wind projected onto the surface of the super- 
giant. The red lines represent the streamlines of the stellar wind 
and the black contour-lines labelled by the corresponding num- 
bers in km - s~ l mark the absolute values of the stream velocity 
in the orbital plane 



cancel each other in our present calculation. However, they may 
be in a different proportion in other systems and with the tem- 
perature structure taken into account. The classical models of 
radial stellar winds in binaries should thus be revised to achieve 
reliable quantitative results. 

The flow of the circumstellar matter in the wider vicinity 
of the mass-losing star can be seen in Figure [4] which repre- 
sents the density, velocity, and streamlines in the orbital plane. 
The colour map on the projection of the tidally distorted disk 
of the star reflects the wind density shown in Figure [2] projected 
from the direction of the rotational axis (its minimum in the polar 
region and the maxima shifted with respect to the region below 
the Lagrangian points can be clearly seen). The colour-map out- 
side the stellar surface represents the density distribution, the ap- 
proximately concentric isocontours show the values of the veloc- 
ity, while the radial curves show the streamlines in the equato- 
rial plane. In general, the density of the stellar wind decreases 
while the velocity increases with the distance from the super- 
giant. The distortion of the streamlines by the Coriolis force is 
obvious and its role in bringing the mass ejected around the L 2 - 
point into the vicinity of the compact companion can be esti- 
mated. (The significance of this effect argued by D. Gies depends 
on the parameters of the system and the regime of the radiative 
acceleration). 

The most spectacular feature in this Figure [4] is the dense 
tail that appears behind the compact companion (the position of 
the companion is marked by the cross). The growth of this tail 
can be seen in Figure [3] Its source is obvious from the stream- 
lines of the wind. The flow of the wind is focused by the in- 
tense gravity field of the black hole, and the mass density in- 
creases in the tail as the streamlines passing the black hole at 
larger distances join it. This effect, which is naturally beyond 
of scope of the models of radial wind, is obviously a version of 



the Bondi-Hoyle-Lyttleton (BHL hereafter) accretion (with dif- 
ferences described in the next Section). In an approximation of 
free-particle motion, this effect is analogous to the gravitational 
focusing of photons and the formation of high-density caustics. 
Our isothermal model facilitates the condensation of the tail, al- 
though a similar increase in the mass density behind the com- 
panion star can be seen in the results of many numerical models 
of mass exchange in binaries. 

An important characteristic of the BHL accretion investi- 
gated already in the original papers of its discoverers but of- 
ten neglected in more recent papers is its instability. There is 
no unique steady- state solution with a strict boundary between 
the region of the tail that falls on the accretor and which escapes 
to infinity. The accumulation of the mass in the tail and the vari- 
ability of the accretion rate with time may play an important role 
in triggering the transitions of X-ray binaries between the radia- 
tively high/soft and low/hard regimes of the accretion. Our com- 
putations also indicate that there are some instabilities in the tail, 
although an increase in the spatial resolution and a sophistication 
of the physical model will be needed to verify or disprove these 
results. 

The close proximity of a luminous X-ray source could influ- 
ence the dynamics of the line-driven wind by producing a change 
in the temperature and ionization structure of the wind material. 
The X-ray heating and photoionization can inhibit the wind ac- 
celeration by reducing population of electron levels available to 
absorb the momentum of the radiation flux from the primary. 
In the case of a very strong X-ray flux, the photoionization ra- 
dius can reach the photosphere of the primary causing the line- 
driven wind to never achieve the escape velocity. It may then 
cease to supply the accretion process with material. Although 
the calculations of the ionization structure were not included 
in the present version of the code, we were able to simulate 
the above-mentioned effect by the assumption that the k param- 
eter from Eq. ([T]) is artificially set to for the wind material, 
which leaves the part of the photosphere facing the X-ray source. 
If the CAK line-driven mechanism is a main source of the repul- 
sive force acting on the wind material, by cancellation of the line 
force, the wind material loses its support and starts falling back 
onto the primary. The situation is captured in Figure [5] The ini- 
tial state (the uppermost left panel) corresponds to the steady- 
state solution depicted in Figure [4] There is an apparent steep 
increase in the density in the region where the line force has 
been artificially switched off. This is caused by the slowdown 
of the wind velocity in this region and also by the higher layers 
of the wind material failing to achieve the escape velocity and 
falling back onto the photosphere. 



4. Bondi-Hoyle-Lyttleton accretion 

The standard theory of BHL accretion assumes that the accreting 
object moves th r ough a homogeneous st atic gas (cf . Hoyle & 
|Lyttleton| ( fT939] ), |Bondi & Hoyle] ( [1944 ), or Edgar] ( |2004) for 
a review). This means that the initial velocities of the incident 
gas elements are parallel, while the stellar wind diverges. This 
difference is essential because the divergence limits the mass rate 
of the accretion into the formed tail (cf. Hadrava & Karas|1984| ). 
Moreover, for an accretor moving in either a perpendicular or 
skewed direction with respect to the source of the wind, the axial 
symmetry of the flow is broken. 

To modify the calculation of BHL accretion in the case of 
the gravitational interaction of a radial stellar wind with an or- 
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is a symmetric matrix describing the dilution of the unperturbed 
wind in the space around the satellite with respect to its rest 
frame. Its projection to the plane perpendicular onto vq 

k k , j w(R) , d (w(R)\ k R l R j j 

A* = P k .a\P\ = —--P* + — — — Pi , (16) 

1 1 J 1 R 6R\ R ) R 1 

where the projector 
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Fig. 5. Evolution of the density distribution in the equato- 
rial plane computed with the radiation hydrodynamic model 
after the line-driven force from the companion-facing hemi- 
sphere was artificially switched off. The initial state (the upper- 
most left panel) corresponds to the steady-state solution depicted 
in Figure [4] The bottom right panel captures the new nearly 
steady-state solution. The white index in the upper right corner 
of the each panel indicates an evolutionary time expressed in 
the units of the orbital periods. The red cross and asterisk show 
the position of the black hole and the centre of the supergiant, 
resp. 



biting satellite, the boundary conditions must be modified. We 
define 



w l (r) = w(r) — 
r 



(12) 



to be the unperturbed velocity of the wind and R l (t) to be 
the track of the satellite. The relative velocity of the wind in 
the instantaneous inertial system of the satellite can be expanded 
in variable x = r - R up to linear terms as 



v l (x) = w l (R + jc) - R l = v l + a l jX J 
where 

^ = w(/?)|-7?' 



(13) 



(14) 



is the velocity of the wind in the rest frame of the satellite col- 
liding with it head-on, and 



, w(R) si d w(R)\R'RJ 



(15) 



P' = 8'- - 

1 1 v 2 



(17) 



yields the apparent divergence of the wind streamlines. For 
a vector z that is perpendicular to the orbital plane (i.e. (zR) = 
(zM) = (zvq) = and consequently Pz = z), we get 



v(R)z k 
R 



(18) 



which means that the dilution of streamlines in this plane cor- 
responds to the proper distance R of their source. On the other 
hand, for a vector y in the orbital plane perpendicular to vq, i.e. 



y = c[R l (R 2 - ^-(RR)) + R l (Rw(R) - (RR))], 
R 

(for which (yvo) = and hence Py = y, but (yR) 

cw(R) ,- P 2 hi 



(19) 



R 



[R Z R Z - (RR) ] * 0), we get 



A\y l 



(20) 



w(R) 



R 2 R 2 - (RR) 



o\2 



d (w(R)\ 



R Rw 2 (R)~2(RR)w(R) + RR 2 dR\ R 



y 



In particular, when the satellite moves perpendicularly to 
the wind (i.e. (RR) = 0) and the wind is not accelerated 
(^ (w(R)) = 0), the apparent divergence of streamlines in the 
orbital plane is reduced as if the distance of their source is in- 
creased by the square of the ratio of the apparent (vq) to the true 
(w) wind velocity 



w 3 (R)y k 
fa 



(21) 



As a consequence, the streamlines in the equatorial plane 
(x,y in Figure [6]) are focused on the x-axis more efficiently than 
the streamlines in the perpendicular plane (x, z). The streamlines 
passing the satellite out of both these planes cross the x, z-plane 
earlier than the orbital plane and are thus not focused exactly 
on the x-axis; the focusing is astigmatic. This conclusion agrees 
with the results of numerical models described in the previous 
Section, which show that the cross-section of the tail is elongated 
in the direction perpendicular to the orbital plane. The more 
complicated structure of the tail will thus need to be treated using 
detailed 3D hydrodynamic modelling to study the time variabil- 
ity of the accretion rate. 



5. Discussion and conclusions 

The increasingly high quality data available owing to improved 
observational techniques at various wavelengths now enable 
the more detailed investigation of the circumstellar matter in 
the interacting binaries. An accurate interpretation of these ob- 
servations requires more sophisticated theoretical models to get 
more precise predictions of the observable quantities and test 
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Fig. 6. Gravitational focusing of the wind in the satellite's 
co-moving frame in the approximation of free-particle motion. 
The thin black lines show asymptotes of the Keplerian hyper- 
bolic trajectories. Owing to an aberration of the initial veloci- 
ties, the trajectories with the same impact-parameter in the or- 
bital plane (x, y) (red in the electronic version) reach the x-axis 
closer to the gravitating body (placed in the coordinate centre) 
than those in (x, z)-plane (blue lines). The particles moving out 
of these planes (violet lines) collide mutually in the (x, z)-plane 
before reaching the x-axis 



observationally the viability of the models and their underlying 
assumptions. 

To yield physically self-consistent and quantitatively reliable 
results, our numerical model of the 3D radiation hydrodynam- 
ics of stellar winds in high-mass X-ray binaries (cf. Negueruela 
2010 for a recent review) presented here needs to be improved in 
terms of several aspects of physics (first of all to model the radia- 
tive transfer of the X-rays from the companion and their interac- 
tion with the wind) and numerical properties (higher resolution). 
However, even in its present form it shows that non-radial mod- 
els are necessary to derive quantitative information about the an- 
gular modulation of the wind and the mass-loss and accretion 
rates with an accuracy of about ten percents. The preliminary re- 
sults of our calculations in which the radiative pressure in lines is 
scaled to account qualitatively for the ionization of the wind by 
X-ray illumination in different states of the compact companion 
show the significant non-linear response of the wind strength to 
the radiative drag. 

The formation of narrow gaseous tails (or 'wakes') behind 
compact objects can be seen in our numerical hydrodynamic 
models of high-mass binaries. Similar results have also been 
found in several previous hydrodynamic simulations of isother- 
mal flows - e.g. |Ishii et aT\ ( |1993| ), while the adiabatic flows 
tend to form wider shock cones - cf. |Matsuda et al.| ( |1991| ). 
(Note that the equations of motion (3) and (2.2), resp., are in- 
correctly typed in these papers.) The physics of the formation 
of this tail is a modification of the BHL- accretion, i.e. a gravita- 
tional focusing of the wind that is basically similar to the gravi- 
tational lensing of photons. Unlike the axially symmetric case of 
the classical BHL-accretion, the stellar wind from a supergiant 
has spatial gradients of both the velocity and the density when 
reaching the orbit of the companion object and its accretion must 
thus be non-axisymmetric. The violation of the symmetry is 
more pronounced if the wind speed does not significantly exceed 
the orbital velocity of the companion. In this asymmetric case, 
the cla ssical approxi mation of the BHL accretion rate is less ac- 
curate ( |Walder|1997| ) and the accretion is unstable (Ruffert 1997] 
19991. There are sophisticated general-relativistic radiation- or 
magneto- hydrodynamic models (cf., e.g., |Zanotti et al.||201 1" 



needed to treat in detail the accretion in the close vicinity of 
the compact component and the generation of the X-ray radi- 
ation. The present models, however, mostly assume the homo- 
geneity of the incident flow. An inhomogeneity in the flow may 
significantly influence the angular momentum of the accreted 
matter, hence both the size of the inner accretion disc and its lu- 
minosity. The inhomogeneity of the wind in the neighbourhood 
of the companion depends on its tidal modulation in the subsonic 
region of the donor star and its non-radial flow in the space be- 
tween the components. These effects have thus to be taken into 
account in self-consistent models of the accretion. 

The existence of these gase ous tails in interact ing binar ies is 
confirmed observationally - cf. D umm et al.| ( |2000] ) or |Dhawan] 
et al.| ( 200 6). The deposition of the mass focus ed into the tail 



and the well-known instability of BHL accretion (Foglizzo et al 
20051 may influence the interplay between the variability of 
the wind mass-loss from the donor star, the hydrodynamics of 
the circumstellar matter, and its accretion onto the compact ob 
ject. It may thus be responsible for different kinds of time vari 
ability for these objects - cf., e.g., Blondin lfe Pope| (2009). 



and references therein) of the BHL-accretion. These models are 
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Appendix A: Model atmospheres of binary 
components 

We now reproduce in a language corrected form the parts of 
Hadrava (1987) relevant to our present study. Figure captions 
missing in the original paper because of lack of space are added 
here. 

The models of stellar atmospheres are based on the assump- 
tion of either their plane-parallel or spherical symmetry. The vi- 
olation of either of these assumptions by tides and rotation in 
close binaries leads to the incompatibility of hydrostatic and ra- 
diative equilibria. An improvement in the model atmospheres in 
this respect is desirable for simulating both the light curves and 
line profile changes. Moreover, the atmospheres of the contact 
components of interacting binaries determine the initial condi- 
tions of dynamics of gaseous streams and in this way influence 
the behaviour of the binary system. 

A.1. Hydrostatic equilibrium 

It is well-known that the stellar configuration in hydrostatic equi- 
librium must be homogeneous across each equipotential. This 
follows from the vector nature of the equation of the hydrostatic 
equilibrium 



pVO + VP = , 



(A.l) 



because the pressure P and density p must be functions only of 
the potential ®, and following the equation of state 



P = P(P,T) 



(A.2) 



the same must be valid for the temperature T = T(<S>). The tem- 
perature distribution is determined by the energy balance, driven 
mainly by radiative transfer. The solution of radiative transfer 
is the main problem in modelling symmetric (i.e. plane-parallel 
or s pheri cal) atmospheres, where the t ange ntial c omponents of 
Eq. ( |A.1| ) are trivial and both Eqs. ( |A.1| ) and ( |A.2| ) can be simply 
integrated. However, if the symmetry is violated e.g. by either 
rotation or tides, the homogeneity imposes an a priori restric- 
tion on the radiation field across different parts of stellar sur- 
face. This restriction cannot generally be satisfied. This can be 
seen from Fig. A.l where the source function S is plotted (full 
line) for the grey atmosphere. If g is increased, the hot inner 
layers are shifted to lower optical depths (dashed line), which 
in agreement with von Zeipel's theorem corresponds to an in- 
crease in the effective temperature (dotted line). On the other 
hand, the surface temperature (for r —> 0) remains constant and 
is thus incompatible with the radiation coming from the interior. 
When the surface layers are heated by absorbing this radiation, 
the hydrostatic equilibrium is violated in the tangential direction 
and tangential radiative transfer takes place. The latter of these 
generally also appears as a consequence of the different radia- 
tion dilution caused by the variable curvature of the equipoten- 
tials. The incompatibility of the hydrostatic and radiative equi- 
libria should give rise to meridional circulations in the atmo- 
sphere. The construction of model atmospheres as a mosaic of 
one-dimensional (especially plane-parallel) models can thus be 
accepted only as a first approximation. Another violation of hy- 
drostatic equilibrium is caused by the existence of stellar wind, 
which is discussed in Section |A.2| Additional complications 
arise when the binary component does not rotate synchronously 
with the orbital motion, particularly when the orbit is eccentric. 
Despite the aforementioned objections, the hydrostatic equilib- 
rium is a valuable approximation, at least for the inner layers 




Fig. A.l. The dependence of source function S in a grey at- 
mosphere on the optical depth r (the full line). In the parts of 
atmosphere with a higher gravity g, layers with a given tem- 
perature are shifted to lower optical depths (the dashed line). 
This corresponds to increase of the radiative flux in large opti- 
cal depths (the dotted line), which gives, however, a temperature 
incompatible at low depths with the homogeneity on equipoten- 
tial. The actual value of emergent flux H for different shifts of S 
in g is shown on Fig.|A.2 




Fig. A.2. The dependence of radiative flux H from an inho- 
mogeneous grey atmosphere in exact hydrostatic equilibrium on 
the local gravity g (in units of its mean unperturbed value). 



of atmospheres. It is thus advantageous to choose a coordinate 
system that fits the corresponding shape of the atmosphere for 
the further investigation. This is achieved using Roche coordi- 
nates, where the first coordinate £i is the effective potential and 
the surfaces £2,3 = const are perpendicular to the equipotential s 
(but generally not mutually perpendicular). The generalisation of 
a Roche potential to a finite mass concentration is also possible, 
but fortunately not essential. 

A.2. Stellar wind 

The density of a gas in thermodynamic equilibrium is an expo- 
nential function of the potential. This is why it decreases to zero 
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and the mass (or optical) depth is finite in upper layers of plane- 
parallel atmosphere, where O — > oo. However, O — > at infinity 
for a spherical (or other finite) configuration and consequently 
p — > Poo > and the mass depth becomes infinite for such a con- 
figuration in equilibrium with its surroundings. If a real star is 
immersed into vacuum, it must evaporate. Eq. ( |A.1| ) of hydro- 
static equilibrium must then be thus replaced by the equation of 
motion 



d t u + (uV)u + 2Q A u + V® + -VP = 

P 

and the equation of continuity 
d t p + V(pu) = . 



(A.3) 



(A.4) 



This system of partial differential equations can be simplified for 
a stationary flow (d t = 0), if the stream lines 



x(s) . 



u = un 



dx 

ds 



(A.5) 



are known (e.g. if they are radial for spherically symmetric prob- 
lem). Eq. \AA\ then becomes 



ds 



ipuD) = 



and its solution is 
h 



(A.6) 



(A.7) 



where h is the flow of gas, which is constant along the stream 
line, and D = D(s) is the cross-section of stream-tube, deter- 
mined by the dilution of the stream lines 



ds 



InD = (Vn) . 



The component of Eq. ( |A.3| ) parallel to n reads 



d u 2 d ^ Id 
— — + — 0+ --P 
as 2 as p as 







(A.8) 



(A.9) 



The density can be substituted here from Eq. ( |A.7| ), and if 
the thermal regime of the flow is known (e.g. where there is 
a polytropic expansion or T = T(s)), the equation for u = u(s) is 
obtained. It can be simply integrated in the form 



2 2 
U k ~ U 



kT , 



kT , u 
+ — In — 

m Uk 



(A. 10) 



= InD + O - O fc = (DeffCv) - cD A 

m 



for an idealised problem of isothermal flow (O^ is an integra- 
tion constant). There are two solutions u(s) < = and 

u(s) > Uk of this equation - see Fig. |A.3| The former is asymp- 
totically hydrostatic, while the latter corresponds to a free fall 
for ® e ff — > -oo. An important feature is that the dilution of 
the stream lines acts as a repulsive force (nozzle) through the first 
term on the r.h.s. of Eq. ( |A.10| ). The above-mentioned solution 
can be matched at the critical point Sk, where O e ff reaches its 
maximum ® e fffe) = and u(st) = u^. Its position can be in- 
fluenced or even dominated by the presence of another repulsive 
force such as radiative pressure or centrifugal force. Thei latter 
is combined with the tidal force in binaries. A critical point thus 
depends on the Roche potential along a particular stream line. 
It determines the velocity in deep layers, where the hydrostatic 




Fig. A.3. Evaporative wind in an effective Roche potential - 
the left an d right upper panels give the left and right sides of 
Eq. ( |A.1Q ), i.e. the effective potential ® e ff as a function of u and 
s, resp. The bottom panel gives the solution u = u(s), the full line 
corresponds to the wind solution, and the dashed line to accre- 
tion. 



equilibrium is approached and the density is the same exponen- 
tial function of potential for each stream line. The flow of mass 
is thus 



h ~ exp 



~~kl 



(A. 11) 



in d iffere nt directions. The results for a radial wind are shown in 
Fig. A.4 for several values of the mass ratio q and ratio 6 = 
of thermal kinetic to gravitational potential energy of a gas par- 
ticle. The Roche potential is thus a potential barrier to the cool 
gas, that can flow only in the vicinity of the L\ -point, although 
the anisotropy of the stellar wind decreases with increasing sur- 
face temperature. In any case, the Roche-lobe overflow and mass 
loss by stellar winds are not independent processes as often as- 
sumed, but only two extreme pictures of the same physical pro- 
cess - an anisotropic stellar wind. 
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■/Q\/CX/C\ 

Fig. A.4. Angular distribution of wind strengths for several val- 
ues of and q. Full lines mark the isocontours of the mean wind 
intensity (dashed lines show values for one and two orders higher 
and lower) in spherical coordinates with poles towards the L\ 
and L 2 points. The latitude y> is drawn in a linear scale on the hor- 
izontal axis; the longitude X measured from the intersection with 
the orbital plane scaled by cos (p increases upward up to the ro- 
tational axis at the top of each drawing. 
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